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ABSTRACT 


Greek symbols 


The RNG-based algebraic turbulence model, with a 
new method of solving the cubic equation and applying 
new length scales, has been introduced. An analysis has 
been made of the RNG length scale which was previously 
reported and the resulting eddy viscosity has been com- 
pared with those from other algebraic turbulence mod- 
els. Subsequently, a new length scale is introduced which 
actually uses the two previous RNG length scales in a 
systematic way to improve the model performance. The 
performance of the present RNG model has been demon- 
strated by simulating the boundary layer flow over a flat 
plate and the flow over an airfoil. 

NOMENCLATURE 

Ad constant (=0.2) 

b constant (=§A<*~^) 

c airfoil chord length 

C wave length cut-off constant (75-200) 

Cj skin friction coefficient 

Cp turbulence model constant (=0.09) 

FkihI KlebanofTs intermittency factor 
H Heavyside function 

l mixing length 

li mixing length from Ref. 2 

/ 2 mixing length from Ref. 4 

Lj integral scale of turbulence 

Re c Reynolds number based on chord length 
Re x local Reynolds number 

Re$ Reynolds number based on local momentum 
thickness 

RNG Renormalization Group 

t time 

TI free stream turbulent intensity 

it stream-wise velocity component 

U mean velocity component in streamwise direc- 
tion 

x coordinate direction along chord line 
y coordinate normal to the wall 
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outer model constant (=0.0068) 
turbulence model constant (= <£*) 


boundary layer thickness 

displacement thickness 

turbulent eddy viscosity 

momentum thickness 

von Karman constant (=0.4) 

kinematic viscosity 

the provisional eddy viscosity 

effective viscosity 

density 

wadi friction parameter (= yj^^) 


shear stress 

Coles pressure gradient parameter 
wave number 


Superscripts 


normalized quantity 
time averaged quantity 
+ law-of-the-wall coordinate 


Subscripts 

0 laminar quantity 

1 inner layer 

o outer layer 

t turbulent quantity 

e quantity at outer edge of boundary layer 
rn^l based on length scale from Ref. 2 

rng2 based on length scale from Ref. 4 

bl based on length scale from Ref. 6 


1. INTRODUCTION 


1 


Studies of the transition from laminar to turbulent 
flow in the boundary layer of an airfoil were performed 
in the 1940’s and earlier, upon observing that extensive 



laminar layers existed on the forward portion of smooth 
wings. However, in recent years only limited numbers 
of turbulence models capable of simulating the laminar- 
turbulent transition have been developed. 

The drag of an object such as a flat plate or an 
airfoil with a smooth surface depends on the location of 
the point of transition where the value of drag suddenly 
changes from low-drag laminar-type flow to high-drag 
turbulent flow. An incorrect estimation of the transition 
point greatly decreases the level of accuracy in predict- 
ing the flow characteristics farther downstream. There- 
fore, the accurate prediction of drag by CFD methods 
requires the correct prediction of the transition location. 

In 1986, Yakhot and Orszag 1 proposed the appli- 
cation of the renormalization group theory to turbulence 
modeling for transition simulation. Later Martinelli and 
Yakhot 2 proposed an algebraic and differential k - e 
model based on the RNG theory of Yakhot and Orszag 1 
and applied the algebraic RNG model to transonic flows. 
Lund 3 investigated the algebraic RNG model and used 
a simple cubic equation for v after imposing reasonable 
assumptions. Kirtley 4 used a quart ic equation based 
on the cubic equation and applied the model to three 
dimensional turbomachinery using Newton’s method to 
find the roots of the quartic equation. 

The main reason for using a quartic equation is 
due to the existence of more than one solution to the 
cubic equation whose locus of real roots consists of two 
branches. By special transformation to the quartic equa- 
tion, a difficulty of finding one non- trivial solution was 
simply eliminated and Newton’s method was used to 
evaluate the solution to the transformed quartic equa- 
tion of Kirtley 4 . Therefore, both branches have been 
retained to avoid a jump discontinuity in the distribu- 
tion of effective viscosity, and the solution was evaluated 
according to some forcing conditions given by Lund 3 
[Figure 1]. Kirtley 4 successfully used a method stem- 
ming from direct evaluation of the cubic equation with 
the jump condition and showed improvements over the 
Baldwin-Lomax model with modification for wake and 
separated flows. 

In the present study, one of the solution branches 
that shows less physically important behavior than the 
other real solution branch has also been eliminated and 
the solution has been evaluated according to a general 
solution formula for a cubic equation which involves 
algebraic manipulation of complex numbers. The re- 
sulting jump discontinuity in normalized viscosity does 
not show any negative contribution to convergence in 
the main flow solver [Figure 2]. The present method is 
compact in size and simple in analysis since it only re- 
quires the evaluation of the cubic equation at each mesh 
point. Due to the straightforwardness of the present 
analysis, the implementation effort for coding has been 


minimized. The main intent of the present study is to 
provide a better understanding of the model’s transition 
behavior and to provide extensive flow results in order to 
check the performance of the model on boundary layer 
transition simulation as well as turbulent simulation. 


2. DERIVATION OF THE MODEL 

The mathematical form of the RNG theory is given 
by Yakhot and Orszag 1 as 


v = i/q 


1 + 


TJ (i A 1.594. 


1/3 


( 1 ) 


where H(x) is a Heavyside function defined as follows: 



for x > 0; 
for x < 0. 


In Equation (1), v and i/q are the effective and kinematic 
viscosity respectively, Ad is a constant (=0.2), and A/ is 
a wave number corresponding to the integral scale of the 
turbulence in the inertial range. A cut-off constant de- 
rived by Yakhot and Orszag 1 is represented by C (=75- 
200), and l is the dissipation rate. 

Lund 3 derived the following cubic RNG formula- 
tion from Equation (1) after a series of simplifications: 



where Lj is an integral length scale, and V\ can be de- 
termined from Prandtl’s mixing length theory: 


Vt 


/ 2 



( 3 ) 


Also by applying the fully turbulent condition, 
that is 


i>t + 1 ~v t for v t > 1 


( 4 ) 
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the integral and mixing length scales ratio can be rep- 
resented as follows: 


* )* = 1 + ■ 11 * 1 


(5) 


v t = l\ 


du 

( *y N 

v 2 du 

dy 

\.l + 0(y/6) y 

) dy 


or 


which leads to the final governing equation for effective 
viscosity not only valid in laminar regime but also in the 
transition/turbulent regime, that is 


i>= [l + H (yi>t - C)] 1/3 
Equation (6) is further simplified as 

i> 3 _ + {C - 1) = 0 

with a positive H function condition, 
vi>f - C > 0 


( 6 ) 


(7) 


( 8 ) 


‘"[i + P(y/sW y 


du 

dy 


(ii) 


Stock and Haase 5 rewrote y §^-| as function of y/6 using 
Coles velocity profile formulation: 


du 

dy 


= ^ [x|<r|n| sin(x|) + <r (12) 


The expression for y | §*j- 1 can be approximated by taking 
its maximum. 


The mixing length for wall bounded flow was taken 
from Martinelli and Yakhot 2 


h = 


ky 


i + P(y/6) 


(9) 


for transition prediction and from Kirtley 4 in turbulent 
flow region. 


l 2 = CpS tanh 


(s) 


( 10 ) 


The reasons for using these two length scales will 
be discussed in Section 3 and 4. 


Q = 7 sin(w^), for 0 < < 1 

0 0 V 


where, the maximum value of <5 can be obtained by 
Newton’s method; 


Q\max — 1-82 at ~ 0.646 (13) 


Using these values, Equation (12) becomes 

U c a 


da 

dy 


■[1 .8211 + 1] 


(14) 


The result of substituting Equation (14) into (11) is: 


3. EVALUATION OF RNG LENGTH SCALE 


_ (0.372) 2 (0,646)6 U e a 

[ 1 + ( 4 . 44 )( 0 . 646 )p * [ L82II + J 


The RNG length scale originating from Martinelli 
& Yakhot 2 [Equation (9)] has been analyzed using 
the Stock k Haase 5 analysis and compared with the 
Baldwin-Lomax 6 model. 

By substituting Equation (9) into Equation (3), 
the provisional eddy viscosity, i.e., the maximum eddy 
viscosity from RNG model, can be written as 


or 


1 C£ 0.006 U e — [1.82H + 1] (15) 

/c 


The same procedure has been performed by Stock 
and Haase 5 for the Baldwin-Lomax 6 turbulence model’s 
outer layer viscosity formulation yielding, 


3 



i/ t>M k 0.0173 U e — [1-82H + 1] (16) 

K 

Using Equations (15) and (16), the ratio of Vt t mg\ and 
VtM is: 

^ZSL = iM. « 0.35 (17) 

i/ t M 0.0173 

From Equation (17), it is shown that the RNG viscosity 
predicted by the length scale of Equation (9) is about 
one third of that from the first Baldwin-Lomax outer 
model of Equation (16) that covers the near wall and 
attached flow regions. Consequently, it can be concluded 
that the length scale of Martinelli and Yakhot 2 predicts 
a much lower value of the eddy viscosity than the length 
scale of Baldwin-Lomax outer model. Also it has been 
shown that this length scale resulted in very low skin 
friction coefficients for flat plate and airfoil flows. 

The length scale of Kirtley 4 [Equation (10)] has 
been introduced in the current RNG model in combi- 
nation of that of Martinelli et al. 2 . The same analysis 
of Stock & Haase 5 can be performed again. Since the 
maximum length scale occurs for |tanh x\ maT = 1 in 
Equation (10), the maximum v% for the model is ob- 
tained from Equation (3). 

*w = (<V)’|f;| ( lg ) 

OT e fS\ 105 

for provisional RNG viscosity. Again, using Equa- 
tion (14) and corresponding y ma x/& value, 

Vt, Tn9 ^Cl6-uAlMU + l] (19) 

y I* 

OT - 

0.0125 VtjfllJML + 1] 

From Equations (16) and (19), 

Ut ’ rn ? 2 - 0 0125 ^ o 72 (20) 

vtju 0.0173 1 ' 

From equation (20), the length scale of Kirtley 4 shows 
a more comparative magnitude of eddy viscosity than 


the previous RNG length scale of Martinelli Yakhot 2 
[Figure 3]. Also, the results of skin friction for a flat 
plate and an airfoil using the length scale of Kirtley 4 
are in better quantitative agreement with experimental 
data than the results obtained with the length scale of 
Martinelli & Yakhot 2 . However, one drawback in us- 
ing the length scale of Kirtley 4 is that it predicts the 
laminar-turbulent transition point too early, especially 
for the flows with relatively low free stream turbulence 
intensities. 

After the test of the RNG model using the length 
scales of Martinelli et al 2 (h) and Kirtley 4 (/ 2 ) for flow 
over a flat plate and over a NACA 0012 airfoil, the pre- 
vious analysis has confirmed that the length scale, /i, 
from Equation (9) underpredicted the skin friction up 
to 30% but it accurately predicted the onset of transi- 
tion from laminar to turbulent. Use of the length scale, 
/ 2 , from Equation (10) provides excellent agreement of 
skin friction in the fully turbulent region but predicts 
the transition point too early. Due to the sudden in- 
crease of the length scale near the wall (0 < y/6 < .2) 
in the boundary layer, the solution sometimes diverges 
just after the model has been turned on. 


4. IMPLEMENTATION OF RNG MODEL 

Based on the analysis of the RNG length scales 
from the previous section, a new length scale which pre- 
dicts the transition point accurately and also produces 
the correct turbulent skin friction coefficient is intro- 
duced. 

Because the transition location (y/6)t r between in- 
ner and overlap layer inside a turbulent boundary layer 
is not known initially and it depends on many factors, 
such as initial free stream turbulence, Reynolds num- 
ber, pressure gradient, surface curvature, and surface 
roughness, the actual generation of a single length scale 
may not be easily determined. Consequently, such an 
attempt was not made in the present implementation. 
However an RNG model is introduced, which makes a 
systematic use of the previous two length scales one at 
a time as follows: 

The length scale, /i from Equation (9), is used 
first to find the accurate transition location and, once 
the positive H function condition is satisfied from Equa- 
tion (8), the length scale, li is replaced by from Equa- 
tion (10) for the region where the RNG has been turned 
on. 

From the solution of Equation (7) and ac comply- 
ing condition Equation (8), the model is not turned on 
until the provisional eddy viscosity reaches the following 
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value: 




*=Vz 


[(C- 1)1 1/3 

2 


( 21 ) 


where the provisional viscosity is calculated from the 
Equation (3) as follows: 


du 
dy 

where 



f Zi, for transition prediction; 

\ / 2 , for turbulent RNG viscosity evaluation. 


where Q is a vector containing the conservation vari- 
ables, 


Q = 


p 

pui 

E 


where Fj vectors represent the invicid flux vectors, 


Fj = 


puj 

pmuj + pSij 

( E + p)tij 


and the viscous flux vectors Gj are 


The above implementation of two length scales can be 
depicted only if one new length scale were used as shown 
in Figure 4 where the new length scale follows the /i for 
y/6 < (V/S)tr and / 2 for y/6 > (y/<5) tr - 

In the actual implementation, unlike the Baldwin- 
Lomax model where the turbulent viscosity is calculated 
for the entire computational domain, the following algo- 
rithm was used in order to avoid unnecessary computa- 
tion for turbulent viscosity in the laminar region: 

Since the transition criterion for v x is known in 
Equation (21), for the laminar region, only the calcula- 
tion of Equation (3) has been necessary and the eval- 
uation of the RNG viscosity through the Equation (7) 
can be bypassed with the solution of v = 1, which is 
the laminar viscosity value. Once the transition has 
been triggered by the condition in Equation (8) and v% 
reaches the critical value in Equation (21), evaluation of 
the RNG viscosity via Equation (7) is performed with 
v% calculated using / 2 . 

With its non-iterative evaluation scheme of the 
RNG viscosity from the cubic equation, [Equation (7)], 
the above length scale implementation shows much im- 
proved computational results relative to the previous 
one using either one of the above mixing lengths. 


5. COMPUTATIONAL RESULTS 

The Reynolds- averaged Navier-Stokes equations 
(continuity, momentum, and energy equations) can be 
expressed in the following conservation law form: 


<9Q 9Fj _ 5Gj 
dt + dxj dxj 


G j = 


0 

T ij 

UkTjk -qj m 


The eddy viscosity concept and the algebraic forms of 
the Renormalization Group turbulence model are used 
for turbulence closure. A simple algebraic eddy viscosity 
model based on the Renormalization Group turbulence 
model is used for closure of the compressible Reynolds 
averaged equations. Two turbulent boundary layer flows 
using the RNG turbulence model with the mixing length 
proposed earlier were solved together with the Navier- 
Stokes equations for compressible flows in full conserva- 
tion form. The numerical method is a finite difference 
discretization with the Beam and Warming 7 approxi- 
mate factorization algorithm to integrate in time to a 
steady state. The calculations were performed on a Cray 
YMP computer at NASA Lewis Research Center. 


5.1 Flow over a flat plate 


Experimental data for the flow over a flat plate 
with a zero pressure gradient can be found in Sohn 
& Reshotko 8 . The Mach number was 0.2 and the 
Reynolds number based on plate length was approxi- 
mately 2.3 x 10 7 with a physical dimension of 16 feet. 
The corresponding computational domain, covering only 
the upper half of the plate, starts from the free stream 
1 foot ahead of the leading edge and extends to the end 
of the plate in the streamwise direction. The domain 
extends 3 feet high from the wall in the normal direc- 
tion. A grid of 111 x 81 points was used with about 
one third of these located within t/ + < 1000. About 15 
grid nodes were put within t/ + < 50 in the normal di- 
rection. The streamwise grid was densely packed near 
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the leading edge for an accurate prediction of the large 
gradients in the vicinity of the wall and the leading edge 
[Figure 5]. 

In order to check the grid-independence and also 
obtain a more accurate prediction of the transition lo- 
cation, the calculation had been repeated with a finer 
grid of 221 x 161 points. No variation of transition lo- 
cation had been obtained. The converged solution was 
obtained after approximately 10,000 iterations. 

As shown in Figure 6, the transition can be lo- 
cated where the non-dimensional boundary layer thick- 
ness along the streamwise direction suddenly increases 
from a constant value. In the figure, the present com- 
putation clearly shows similar trend and agrees quali- 
tatively with experimental data of Hanson 9 . The com- 
putational results for the skin friction coefficient and 
the shape factor are compared with experimental data. 
Figure 7 shows the local skin friction coefficient and Fig- 
ure 8 shows the shape factor distribution. The solid line 
represents the computational results using the present 
RNG model, and the symbols represent the experimen- 
tal data under different free stream turbulence intensi- 
ties. The analytic results for laminar flow and a semi- 
empirical correlation for turbulent flow are also plotted. 
The results are in good agreement with the experiment, 
and it is shown that the present RNG model, without 
much dependence on empiricism, can accurately mimic 
the detailed phenomena of laminar-turbulent transition. 
Figure 9 shows the turbulent boundary layer velocity 
profile neax the end of the plate showing excellent agree- 
ment with a correlation of Muster 10 in both inner and 
outer layers. 

Setting up a convergence criteria for RNG simula- 
tion prior to computation is not straightforward since, 
as plotted in Figure 10, the final residual at convergence 
is not the minimum value throughout the computation 
due to a residual jump after the RNG model has turned 
on. In the figure, two curves are depicted from two sim- 
ulations using different time step control; namely, time 
accurate simulation with constant time-step size (At=5) 
and local time stepping simulation with variable time- 
step size (maximum time-step size allowed during the 
run, A imax, is set to 20). The same skin friction results 
from the two cases were obtained after 15000 iterations 
for At =5 and 6000 iterations for At max =20. The ir- 
regular convergence pattern after the model turned on 
for At mar =20 is due to larger time step and/or the lo- 
cal time stepping function of current code which also 
started after RNG has turned on. Meanwhile, the time 
accurate execution shows more regular pattern. In both 
cases, the sudden jumps in the residual are due to the 
jump discontinuity of eddy viscosity in RNG formula- 
tion. It should be noted, however, that the local distur- 
bance due to switching on the model at transition region 


does not deteriorate other regions and overall computa- 
tion remains stable. 

5.2 Control of transition 

The present RNG model uses the analytic rela- 
tion of the S/% tmax in its length scale when determin- 
ing the provisional eddy viscosity. It has been observed 
that, by increasing or decreasing 6 /y ma x from its stan- 
dard value (=1.548), the change of the transition loca- 
tion can be simulated numerically as if the flow is in a 
different environment which affects the transition. The 
two computational results plotted in Figure 11 and 12 
as dashed and dotted lines are RNG results obtained 
by using different boundary layer thickness parameters, 
6/ y m ax- The results obtained by a 25% increase and de- 
crease of SI y m ax show good qualitative agreement with 
experimental data obtained with higher and lower free 
stream turbulence intensities, respectively. This con- 
firms that, as mentioned earlier, this parameter can be 
used to predict the transition location under different 
flow environments which affect the transition location. 

5.3 Flow over a NACA 0012 airfoil 

Another calculation was performed for the flow 
over a NACA 0012 airfoil with zero lift in order to test 
the present turbulence model for a flow with a pressure 
gradient. The experimental data considered here can 
be found in Becker 11 . The air speed was 230 mph and 
the Reynolds number based on the airfoil chord was ap- 
proximately 1.0 x 10 7 with a chord length of 5 feet and 
a maximum airfoil thickness of 12% chord. The corre- 
sponding computational domain covers the entire airfoil 
starting from the free stream at one chord length ahead 
of the leading edge and extending to 5 chord lengths 
in the streamwise direction. The domain extends one 
chord length from the wall in the cross-stream direc- 
tion. A grid consisting of 358 x 100 points was used, 
with about one- third of these located within y + < 1000 
in cross-stream grid nodes. The streamwise grid has 
also been densely packed near the leading edge for an 
accurate prediction of the large gradients in that region 
[Figure 13]. 

The grid-independence was checked by refining the 
grid in both directions. Two hundred more grid points 
were added in the streamwise direction. In the direc- 
tion normal to the stream, more grid nodes have been 
included near the wall which reduces the value to 
below one. The converged solution was obtained after 
approximately 6,000 iterations. 

The computational results for the skin friction co- 
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efficient and boundary layer velocity profiles were com- 
pared with experimental data. Figure 14 shows the 
local skin friction coefficient, and Figure 15a through 
Figure 15c show the boundary layer velocity profiles at 
different streamwise locations. The solid line represents 
the computational results using the present RNG model, 
the symbols represent the experimental data. Since the 
experiment was performed in a high speed wind tun- 
nel with very low free stream turbulence (~ 0.01%), 
the use of the standard value of S/y ma x [Equation (13)] 
resulted in the transition location being predicted at 
4% chord while the experimental value is approximately 
15% chord. This is consistent with the fact that, from 
flat plate results, the standard value of 6/y ma z predicts 
the experimental transition location with relatively high 
turbulence intensity which usually triggers early tran- 
sition (see Schlichting 12 ). The transition location has 
been delayed by decreasing the 6/ymax value following 
the earlier analysis, and the transition location of 15% 
chord was obtained with 6/y ma x 40% lower than the 
standard value. 

In Figure 16, the pressure coefficient has been com- 
pared with experimental data for the upper surface of 
the airfoil. The results indicate good agreement with 
the experimental static pressure distribution along the 
streamwise surface. 


6. CONCLUSIONS 

The renormalization group algebraic turbulence 
model was introduced in order to test its performance on 
boundary layer transition simulation. The new imple- 
mentation method enables the use of two different length 
scales in systematic way in order to predict transition lo- 
cation better and also produce accurate turbulent flow 
quantities. The followings are detailed conclusions of 
the present analysis: 

1) On the basis of the boundary-layer results, it 
was shown that dropping one solution branch of the cu- 
bic equation derived from the RNG theory does not af- 
fect the performance of the model nor the convergence 
of the Navier-Stokes calculation, 

2) The analytic study and computational results 
show that the RNG length scale introduced by Mar- 
tinelli and Yakhot 2 provides excellent predictions of the 
transition location, but it underpredicts the magnitude 
of eddy viscosity resulting in a lower skin friction coef- 
ficient. 

3) The RNG length scale introduced by Kirtley 4 , 
provided excellent predictions of fully turbulent quanti- 
ties, but indicated an early transition location and rapid 


transition to fully turbulent flow which can cause numer- 
ical instability. 

4) A length scale which uses the two previous RNG 
length scales in a systematic way to preserve the merits 
of the two scales is introduced. This new length scale 
improved the model’s performance in predicting transi- 
tion as well as the turbulent quantities as shown through 
a flat plate and an airfoil calculation. 

5) For separated flows or strong adverse pressure 
gradient flows, the analytic study of the RNG length 
scale shows the need to modify either the RNG length 
scales, or the method to evaluate the provisional eddy 
viscosity, in order to avoid too large of an eddy viscosity 
value. 

6) The present model’s capability of predicting 
the transition location under different flow conditions 
has been introduced and delivers satisfactory results for 
flows with different free stream turbulence intensities. 

7) The fact that the present RNG model does con- 
tain less experimental constants or adjustments, and 
yet gives accurate transition predictions as well as fully 
turbulent predictions with its simple algebraic formu- 
lation makes the model general, realistic, cost-efficient 
and comparable to the other low-Reynolds number two 
equation turbulence closure models. 
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APPENDIX 
Implementation to PARC code 

The renormalization group theory (RNG) turbu- 
lence model in the present paper has been successfully 
implemented into the two- and three-dimensional PARC 
code with an additional subroutine RGTURBW and one 
function subroutine EDR. The user can select the RNG 
turbulence model by setting the namelist input param- 
eter IMUTUR=4. Also in order to correctly predict the 
transition location under different flow conditions, such 
as free stream turbulence or surface roughness etc., the 
user can set a new namelist input parameter D20YMX 
within ± 40% of its default value (=1.548) to move the 
transition point. The default value of D20YMX can sim- 
ulate a turbulence intensity of about 2% for the flow over 
a flat plate with zero pressure gradient. The user can 
set a lower value than 1.548 for D20YMX to simulate 
the flow with a lower free stream turbulence. 
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Figure 1. RNG viscosity, v as function of provisional viscosity, v t 
via cubic equation for C=75. 



Figure 2. Modified RNG viscosity, v as function of provisional viscosity, v t 
via cubic equation for C=75. 
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Figure 3. Comparison of RNG viscositys with outer-eddy viscosities of Baldwin-Lomax turbulence model. 



Figure 4. Length scales used in the present model and systematic transition from 
one length scale to another inside a turbulent boundary layer. 
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Figure 15a. Boundary layer velocity profiles in laminar region along with 
experimental data. - NACA0012 Airfoil 



Figure 15b. Boundary layer velocity profiles in transition region along with 
experimental data. - NACA0012 Airfoil 
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Figure 15c. Boundary layer velocity profiles in turbulent region along with 
experimental data. - NACA0012 Airfoil 
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Figure 16. Computed Static-pressure coefficient along with experimental data. 
-NACA0012 Airfoil 
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